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Pablo LagungQ 

Department of Astronomy and Astrophysics 
Penn State University 
University Park, PA 16802 

Loop quantum cosmology homogeneous models with a massless scalar field show that the big-bang 
singularity can be replaced by a big quantum bounce. To gain further insight on the nature of this 
bounce, we study the semi-discrete loop quantum gravity Hamiltonian constraint equation from the 
point of view of numerical analysis. For illustration purposes, we establish a numerical analogy 
between the quantum bounces and reflections in finite difference discretizations of wave equations 
triggered by the use of nonuniform grids or, equivalently, reflections found when solving numerically 
wave equations with varying coefficients. We show that the bounce is closely related to the method 
for the temporal update of the system and demonstrate that explicit time-updates in general yield 
bounces. Finally, we present an example of an implicit time-update devoid of bounces and show 
back-in-time, deterministic evolutions that reach and partially jump over the big-bang singularity. 
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Loop Quantum Cosmology (LQC) provides a frame- 
work to explore the dynamics of the Universe deep inside 
the Plank regime [2: 13 ■ ^ n LQC, quantum evolution is 
no longer dictated by the second order Wheeler-DeWitt 
(WDW) differential equation, but instead the dynamics 
is governed by a second order finite difference equation, 
the so-called LQC Hamiltonian constraint. The WDW 
and LQC equations differ only when quantum geome- 
try effects cannot be ignored. A focus of LQC has been 
the problem of the resolution of the big-bag singularity. 
Early studies [1, H, Q showed that the LQC difference 
equation remains deterministic and well defined at the 
classical singularity, thus allowing evolution that jumps 
over the big-bang. Recent studies @, H, H[ of spatially 
flat, isotropic models with a massless scalar field (j> have 
also shown that, under the LQC framework, the big-bang 
singularity is replaced by a big quantum bounce. 

For these models, the LQC Hamiltonian constraint can 
be written as @ 

d 2 ^{v) = \kB-\v){ 

C(v + 2v ) [9(v + 4v ) - *(«)] - 
C(y-2v )[*(v)-9(v-4o )]}, (1) 

where B(v) — 3 3 |v| | \v + Vo] 1 ^ 3 — \v — Vo] 1 ^ 3 ^ is related 
to the eigenvalues of the inverse volume operator [l(| and 
C(v) — \v\ \\v + v a \ — \v — v \\ /4. Units are such that 
k = 8ttG. In Eq. (TT]) , it is implied that \& also depends on 
<t>. The discrete independent variable v is the eigenvalue 
of the volume operator, with the discretization scale v Q 
fixed by the area gap. The scalar field <f> is used as an 
internal time. Starting with a state that is semi-classical 
at late times, evolutions backwards in time using Eq. ([1]) 
yield a turning point, a quantum bounce, before reach- 
ing the classical singularity at v — 9]. The turning 
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point takes place when the total energy density reaches 
a critical value p crlt < p P i an ck- 

The goal of this Letter \s to gain further insight, from a 
numerical analysis point of view, about the nature of this 
quantum bounce. We will show that the observed turning 
point is strongly tied to the method to carry out the 
temporal updating of the system. Specifically, explicit 
time-update of the semi-discrete Eq. ([1} , such as the 4th- 
order Runge-Kutta updating used in Refs. @, H, El , will 
yield a bounce. One can draw a numerical analogy in 
which the bounce can be understood in the context of 
spurious reflections from finite difference discretizations 
of wave equations in nonuniform grids [ll[ or equivalently 
from numerical solutions to wave equations with varying 
coefficients. Bounces due to discreteness have been also 
obtained in the study of Hawking radiation on a falling 
lattice [iH]. 

Motivated by understanding the differences between 
various ambiguities in the constraints for the full theory, 
we will also show that the bounce could be avoided by 
introducing an implicit time-update, which is formally 
equivalent to adding ad hoc higher order terms to the 
LQC Hamiltonian constraint. With this implicit update, 
it is possible to evolve backwards in time a state that 
reaches the classical big-bang singularity. At that point, 
the wavefunction ^ gets partially reflected and partially 
transmitted through the singularity. The resulting evo- 
lution is deterministic and thus also resolves the problem 
of the big-bang singularity. 

We first introduce the standard mesh-index notation 
used in finite differences, namely Vj = j Av with Av = 
Av . In this notation, the semi-discrete Eq. ([T]) takes the 
form 

tifij = ^B~ X D V {CD V ^ ) (2) 

where B = (Slvl^D^v] 1 / 3 ) 3 /^ and C = \v\D v \v\ with 

D v fj = {f j+i/2 - /j-i/2)/Ai> and D v fj = (/ 3+ i/ 4 - 
fj_iu)/(Av/2) discrete finite difference operators. In 
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the limit Av — > 0, both of these operators become d v as 
well as B — > and C — > Also in that limit, Eq. @ 
becomes d; 2 ,'!/ = 3nv/2d v (vd v '&); that is, one recovers the 
WDW equation in the v oc a 3 representation where a is 
the scale factor. 

To simplify our analysis, we will work with a version 
of Eq. © in which B w l/bj| and C « |u|; that is, 
9 2 #j = (3K|w|/2)D t ,(|w|D„* J ). The differences between 
this equation and Eq. @ are of 0(Av 2 ) and do not af- 
fect our conclusions regarding the big quantum bounce. 
In addition, since ultimately the mean trajectories of 
states in the <f> — v plane are given by the characteristic of 
= (3K\v\/2D v )(\v\D v ^fj), we will concentrate our 
attention on its principal part, namely = c 2 D^^ 

where c = y/3n/2\v\. 

The starting point of our analysis is applying a time- 
Fourier transformation ^j((f>) — e~ lu "^"^ 3 (uj)duj. 
Thus, — oj 2 ^ j = c 2 D 2 ^ j. Next, we insert the plane wave 
solutions tyj±i = e ±tfcA "\['j with wavenumber k. For 
Avuj/2c < 1, this yields the following dispersion rela- 
tionship 



Av u> 
2c 



± sin 



Av k 



(3) 



For Avlu/2c > 1, the wave is exponentially damped. 
From this dispersion relation, the group velocity V g = 
duj / dk is given by 



Vg = ±c cos 



Aw k 



= ±c 



AVU 

2c 



1/2 



(4) 



It is then clear that the states, solutions of the semi- 
discrete equation, will have different dynamics from 
those of the continuum WDW equation. WDW states 
have a group velocity ±c and follow characteristics v = 
w ^ e ±Y / 3K72(0-0») w jtjj v ^ anc j integration constants. 
On the other hand, LQC states have group velocity (|U) 
\V g \ < c and mean trajectories v = ^/2[e ± V 3 ^ (<#, " <#, * ) + 



with £ 



and 



Vb = Avuj/\/Qk. The value Vb is the location where the 
group velocity vanishes. This is where the characteristics 
reverse direction, the big quantum bounce. The condi- 
tion (Avuj/2 c) 2 = 1 is precisely the condition pj p CT n = 1 
for the quantum bounce derived in Ref. [9( after identi- 
fying lu with P^, the momentum conjugate to "time." 

Applying an identical analysis to the version of LQC 
Hamiltonian constraint in Ref. @, HI, one arrives to 
Eq. (U) with v replaced by the triad p and setting 
c = y/2n/3\p\. The location of the bounce is now at 
Pb — Ap u> J \/8kJ3. Since v oc p 3 / 2 oc a 3 , when translated 
to the a representation, the turning points Vb and pb oc- 
cur at different values of a, or equivalcntly, for different 
total energy densities p cr u- This "coordinate transfor- 
mation" was the key ingredient in Ref. @ for having the 
bounce occurring at p crit < ppianck- 



It is very important to emphasize that the derivation of 
V g and its turning point condition did not depend on any 
specific time-update (e.g. Staggered Leap-Frog, Crank- 
Nicholson, Runge-Kutta, etc). The only two assumptions 
used were the "spatial" discretization in Eq. (TT]) or ([2]), 
i.e. centered second order accurate finite differences, and 
that the time-update is explicit; that is, the numerical 
approximation to d^Sj was assumed to depend only on 
j. In other words, the Lh.s. of Eq. (fTJ) or J2]) only de- 
pends on Vj . From the LQC point of view, the explicit 
time-update arises because the inverse scale factor oper- 
ator used in deriving Eq. (TT]) is diagonal. Also important 
is to keep in mind that, since the truncation errors from 
approximating B m l/\v\ and C « \v\ are small and only 
become relevant in the immediate vicinity of the classi- 
cal singularity, the bounce observed from the LQC equa- 
tion would be identical to that from solving numerically 
the WDW equation with second order center difference 
approximations and explicit time-integration. In this re- 
gard, what distinguishes LQC and WDW equations is the 
specific forms of the spatial operators, continuum opera- 
tors for WDW and discrete for LQC. 

A possible way of avoiding the bounce in the discrete 
WDW equation is to perform the coordinate transforma- 
tion a oc lna oc lnv 1 / 3 . With this transformation, the 
WDW equation becomes = c 2 9 2 \l/ with c = 
Using a uniform mesh in a, the group velocity in this case 
is identical to ((4]) with Av replaced by Aa. Thus, pro- 
vided the initial state satisfies Aaw/2c < 1, the group 
velocity will not vanish. The same coordinate transfor- 
mation can be applied to the LQC Eq. (fTJ); however, be- 
cause the LQC theory requires keeping At; constant, the 
mesh in the a coordinate will not be uniform. It is in 
this context that one can view the bounce as spurious 
reflections due to nonuniform grids [Til ]. 

So far, we have only revisited the semi-discrete LQC 
Eq. (TT]) under the numerical analysis eyepiece. In partic- 
ular, we have demonstrated that the nature of the quan- 
tum bounce is not entirely due to the specific form of the 
"spatial" operator in the r.h.s. of the equation, which is 
the main difference between the LQC and WDW equa- 
tions. Also very important, and closely intertwined for 
the appearance of the bounce, is the explicit time-update 
implied by Lh.s. of the difference equation. If time- 
updates play a crucial role, one could ask whether other 
time-updates could radically change the dynamics, and 
in particular the bounce, while keeping fixed the spatial 
discretization. The answer is affirmative. 

Let us consider the following modification to Eq. 



M*8$9j = —KB~ l D v {CD v ^j) 



(5) 



where M v fj = (/j+1/2 + /j-i/2)/2- This is an example of 
an implicit time update belonging to the class of Keller- 
Preissman 602; schemes [Tl|. When expanded, the Lh.s. 
of the equations reads (d^ j+1 + 2<9|\Pj + 9?\Pj_i)/4. 
Which means that one cannot update the value of 
without the updated values of ^j±i- Implicit time- 
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updates could arise if one considers non-diagonal terms of 
the inverse volume operator 0, [HI, [Hj] • Although these 
non-diagonal terms are likely to only become relevant 
near the classical singularity. 

An interesting observation is that Eq. ([5]) can be 
rewritten as 

d^^^B^D^CD^^+^j D 2 v ^j. (6) 

In the semi-classical regime, where 0(Av 2 ) terms can 
be ignored, Eq. ^ exhibits the same explicit temporal 
update as the standard LQC equation. On the other 
hand, near the classical singularity, one could apply D 2 
to Eq. © and use the resulting equation to eliminated, 
recursively, in the r.h.s of Eq. the terms with d^Sj in 
favor of higher order, spatial finite difference terms. How- 
ever, in order to arrive to an equation that completely 
avoids the bounce, one would have to include enough 
higher order terms, that are effectively equivalent to an 
implicit update. 

To derive the dispersion relation and group velocity 
for the modified LQC difference equation, we concentrate 
once again on the principal of Eq. ([5]) and perform a time- 
Fourier transform. One obtains —lo 2 M 2 ^ j = c 2 D 2 ^>j 
with c = ^/Zk/2\v\. After substitution of the plane-wave 
solution, the dispersion relation in this case takes the 
form 



Aw lu ( Av k 

= ± tan 

2c V 2 

From which the group velocity reads, 
' Av k 



V g = ±c sec 2 



±c 



f 



Av uj 
2c 



(7) 



(8) 



Clearly this type of time-update does not have a bounce 
or turning point. The group velocity does not vanish 
for \v\ =^ 0. The characteristics in this case are v = 



/6S(0-0„) 211/2 



At v = 0, the group velocity 



becomes infinite because c = 0. 

Fig. [U shows the dispersion relation at an arbitrary 
point 1 17 j ^ for the WDW equation (solid line), the LQC 
equation with explicit time-update (long dashed line) and 
the modified LQC equation with implicit time-update 
(short dashed line). It is evident that for explicit time- 
updates, there always exists a frequency for which the 
group velocity V g = dui jdk vanishes and triggers a mode 
with the opposite group velocity. The reason for the total 
reflection can also be understood in terms of the left- and 
right-moving fundamental modes. These modes couple 
because Eq. ([5]) is a difference equation with non-constant 
coefficients [11]. The coupling also occurs, even in the 
case of constant coefficients, if the mesh is non-uniform 
or discontinuous [HI- The coupling is such that a pure 
mode will be completely reflected at the point where V g 
vanishes [l5[. Also evident in Fig. [T] is that implicit time- 
updates of the type in Eq. do not trigger reflections 
because their dispersion relation is monotonic. 



1.5 




kAv/2 



FIG. 1: Dispersion relation at an arbitrary point \v\ 7^ 
for the WDW equation (solid line), the LQC equation with 
explicit time-update (long dashed line) and the modified LQC 
equation with implicit time-update (short dashed line). 



Fig. [3 depicts an example of the characteristics of the 
continuum WDW equation (solid line), the LQC equa- 
tion with explicit time-update (long dashed line) and the 
modified LQC equation with implicit time-update (short 
dashed line). The constants of integration were chosen 
so the three characteristics have the same starting point 
when evolved backwards in time. 

Fig. [3] shows the evolution of a coherent state so- 
lution to Eq. ([5]) using a mid-point trapezoidal time- 
update. The initial data at <j>* = is centered at 
v* = 2 x 10 4 ti o , has momentum P^ = 10 4 and is 
given by * = exp [-(Z/er) 2 + i P^ Z] where Z = <p - 
</>* ± ^2/3 Kin (v/v*) and a = ^2/3 kAP^/P^ with 
AP^/P^ — 7.5 x 10~ 2 . With these parameters, if an 
explicit time-update were to be used, the bounce would 
take place at v^ = 3257 v a . 

At early times when » > dj, the mean trajectory 
of the state follows the semi-classical trajectory v = 
v ^ e y/3K/2{4>-4>,) un tii v jg comparable to v^. At this point, 
instead of approaching a bounce, the state evolves to- 
wards the singularity, reaching it in a finite time. As ob- 
served from Fig.[3j when the wavefunction reaches v = 0, 
it gets partially reflected and partially transmitted across 
the classical singularity. This partial reflection and trans- 
mission is due to the coupling between left- and right- 
moving modes at v — [151 ]. 

Applying tools commonly used in numerical analysis, 
we have investigated those aspect in the LQC finite differ- 
ence equation related to the bounce observed in homoge- 
neous models of massless scalar fields. We find that, be- 
cause the semi-discrete LQC equation has non-constant 
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FIG. 2: An example of a characteristics of the continuum 
WDW equation (solid line), the LQC equation with explicit 
time-update (long dashed line) and the LQC equation with 
implicit time-update (short dashed line). Constants of integra- 
tion were chosen so the three characteristics have the same 
starting point when evolved backwards in time 




FIG. 3: Evolution backwards in time of a coherent state using 
an implicit time update. The initial data at 0* = is centered 
at v* = 2 x 10 4 V and has momentum = 10 4 . 

coefficients, explicit time-updates will typically couple 
the left- and right-moving fundamental modes and trigger 



reflections at the point where the group velocity vanishes. 
Explicit time integrations in LQC are tied to having a di- 
agonal operator in the gravitational part of the matter 
Hamiltonian for models in which one uses a single matter 
field as internal time [1 61 ] - In addition, we have investi- 
gated ad hoc modifications to the standard LQC finite 
difference equation that avoid the big quantum bounce 
while preserving a deterministic evolution across the big- 
bang singularity. The changes focused on replacing the 
explicit time update with an implicit scheme. Our nu- 
merical experiments showed that a semi-classical state at 
late times can be evolved backwards in time and reach 
the classical singularity in a finite time. At that point, 
the state gets partially transmitted. An interesting im- 
plication of this result is the outcome from an evolution 
forward in time of a right-moving (increasing v direc- 
tion) state initially at the "other side" of the classical 
singularity. When the state reaches the classical singu- 
larity, it will yield a reflected wavefunction that remains 
on the other side of the singularity and a transmitted 
wavefunction emerging across the classical big bang into 
the physical sector. 

Our analysis and numerical experiments were aimed 
at investigating the features in LQC homogeneous mod- 
els with a massless scalar field that are related to the 
occurrence of a bounce. A complementary investiga- 
tion, based on effective perturbation theory around a free 
scalar model, can be found in p7| . From a general per- 
spective, our findings agree for dynamics which allows 
wave packets to reach small scales. Additional studies are 
needed to investigate whether more complicated models 
involving the issue of time could introduce modification 
similar to those considered here, in particular instances in 
which the LQC difference equation includes non-diagonal 
elements of the inverse volume operator. Those terms 
could have the potential of playing an important role 
near the big-bang singularity, significantly modifying or 
even preventing the big quantum bounce. 
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